
% init_data.m
global bcL bcR hx hx1

%----------------初始条件
xa = 0;
xb = 2;
u0 = @(x) sin(x*pi);
bcL = 1;
bcR = 1;
tend = 0.8;
M=200;
%----------------
%[xa,xb]是求解域
%u0是初值
%tend是终止时间
%bcL和bcR是边界条件
%hx是单元的长度
%hx1是单元长度的一半方便计算guass积分
%ureal是用来储存初值(真解)在每个积分点的值
%对于固定的i,有ureal(i,:),真解就是i个单元上所有积分值
%----------------------------------------------
hx = (xb - xa)/Nx;
hx1 = 0.5*hx;

Xc = xa + hx1:hx:xb - hx1;

ureal = zeros(Nx,NumGLP);
for i = 1:Nx
    for j = 1:NumGLP
        ureal(i,j) = u0(Xc(i) + hx1*lambda(j));
        %计算初值在每个积分点上的值
    end
end